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The logarithmic velocity profile is the most important statistical law of classical turbulence affected 
by channel walls. This paper demonstrates numerically that the logarithmic velocity profile of a 
superfluid flow appears in quantum turbulence under pure normal flow in a channel. We investigated 
the configuration and dynamics of an inhomogeneous vortex tangle affected by the walls, and found 
the characteristic behavior of the log-law. 
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The most important statistical laws in the field of clas¬ 
sical turbulence are the Kolmogorov -5/3 law of bulk 
turbulence and the log-law of channel wall [l|. The Kol¬ 
mogorov -5/3 law of energy spectrum has been confirmed 
numerically 0-0 and experimentally [E| in quantum tur¬ 
bulence (QT). This paper is the first report of the log-law 
in QT. 


QT essentially consists of quantized vortices that are 
stable topological defects arising from quantum conden¬ 
sation. QT is currently one of the most important prob¬ 
lems in low temperature physics, and is studied inten¬ 
sively in superfluid helium and atomic Bose-Einstein con¬ 
densates [a, 17| . Almost all studies on QT have been de¬ 
voted to its bulk behavior. This study numerically in¬ 
vestigated the boundary effect of QT and confirmed the 
log-law, namely the logarithmic velocity profile, which is 
established in the field of classical turbulence 181 . 


A traditional system of QT is thermal counterflow in 
superfluid 4 He Thermal counterflow is character¬ 

istic of the two-fluid model that describes the system as 
a mixture of viscous normal fluid and inviscid superfluid. 
The relative motion of two fluids is driven by heat injec¬ 
tion in a channel. When their relative velocity exceeds a 
critical value, a tangle of quantized vortices appears and 
grows to make the superfluid flow turbulent. This sce¬ 
nario was confirmed experimentally 1^, U| and numeri¬ 
cally 12j,[13]. However, most numerical studies assumed 
that the flow profile of the normal fluid was uniformly 
laminar because no useful information was available on 
what actually happens to the normal fluid in the channel. 

Recent visualization experiments have changed the sit¬ 
uation by observing an unexpected inhomogeneous pro¬ 
file of thermal counterflow. By using a laser-induced flu¬ 
orescence technique with metastable He]] molecules, Guo 
et al. 141 and Marakov et al. [HI observed an inhomo¬ 
geneous velocity profile of a normal fluid component in a 
square channel. Motivated by this experiment, some re¬ 
searchers numerically studied the inhomogeneous vortex 
tangle between two parallel plates under the prescribed 
Hagen-Poiseuille flow 16 h 19|. and in a square channel 
under the prescribed Hagen-Poiseuille and tail-flattened 


normal fluid flow J20|. Vortices near the walls were denser 
than those in the central region, forming a superfluid 
boundary layer. Yui and Tsubota calculated the flow 
profile of a superfluid due to the inhomogeneous vortex 
tangle. 

This result reminds us of the logarithmic velocity pro¬ 
file of turbulence in a classical viscous fluid. Here, it 
would be useful to briefly describe the well-known log- 
law in classical fluid dynamics [lj. Let us consider a 
turbulent flow between two parallel plates. We take the 
direction of the flow as the a:-axis and the plane of the 
wall as the a,’z-plane so that y is the distance from one 
wall. The log-law states that the mean value u of the x 
component of the turbulent velocity obeys 


u = —(logy + const.), 
n 


(1) 


where k ~ 0.4 is known as the Karman constant and 
v* is a characteristic velocity. This log-law is obtained 
from a physical picture of the boundary effect when the 
viscosity is negligible. The fluid momentum is transferred 
from the central part of the channel toward the wall as a 
constant flux per unit time. The transferred momentum 
is dissipated by the fluid viscosity near the wall. The 
log-law has been confirmed in classical turbulence in a 
channel 0. 


In this paper, we show that the mean velocity of a 
turbulent superfluid in a channel obeys the log-law. In 
contrast to classical cases, the statistical behavior should 
be reduced to the configuration and dynamics of quan¬ 
tized vortices. To elucidate the effect of the channel walls 
more clearly than our previous study [HI , we considered 
QT in a pure normal flow between two parallel plates 

In a vortex filament model, a quantized vortex is rep¬ 
resented by a filament passing through a fluid and has a 
definite vorticity 21]. This approximation is very suit¬ 
able for He II because the core size of a quantized vor¬ 
tex is much smaller than any other characteristic length 
scale. At the zero temperature, the vortex filament moves 
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with the superfluid velocity 

v s = v SjU1 + v Sjh + t) Sia , ( 2 ) 


where v SiUJ is the velocity field produced by vortex fil¬ 
aments, t> Si b is that produced by solid boundaries, and 
v sa is the applied flow of the superfluid. The filament is 
represented in parametric form as s = s(£,f), where £ is 
the arc length along the filament. The velocity field 
is given by the Biot-Savart law: 


(r) = 


7 f (si — r) x ds i 


4-7T 


Si - r[ 


(3) 


where 7 = 1.00 x 10 - ' m 2 /s is the quantum of circula¬ 
tion and the integration is performed along the filament. 
This paper addresses the full Biot-Savart integral fl3 [. 
The velocity field v s ^ is obtained by a simple procedure; 
it is the field produced by an image vortex that is con¬ 
structed by reflecting the filament onto the surface and 
reversing its direction. The dynamics at finite tempera¬ 
tures includes the mutual friction; the velocity of a point 
s on the filament is given by 


s = v s + as' x (v n - v s ) - a's' x [s' x ( v n - u s )], (4) 


where a and a' are the temperature-dependent coeffi¬ 
cients and s' = ds/dt; is a unit vector along the filament. 



FIG. 1. (Color online) (a) Schematics of simulation of pure 
normal flow. The normal fluid component flows along the 
x-axis with a Poiseuille profile between two parallel plates, 
whereas the superfluid component has no external flow. The 
periodic boundary condition is applied to the x and z direc¬ 
tions and solid boundary condition to the y direction, (b) 
Vortex line density as a function of t under pure normal flow 
with tJ n = 0.9 cm/s. QT develops to a statistically steady 
state. 


In this paper, we address pure normal flow 0 - As 
shown in Fig. Oja), the normal fluid component flows 
between two parallel plates, whereas the superfluid com¬ 
ponent has no external flow, i.e., x Sia = 0. The flow 
direction is along the x-axis. The solid boundaries are 
applied at y/D = 0 and 2 with the half-width D of the 
channel. The normal flow is prescribed to be a Poiseuille 
profile: 


v n = u 0 


1 - 


V - D 

D 


(5) 


where uq is a normalization factor and x is a unit vector 
in the x direction. 

Simulations were performed under the following con¬ 
ditions. We discretized the vortex lines into a num¬ 
ber of points held at a minimum space resolution of 
= 8.0 x 10~ 4 cm. Integration in time was achieved 
using a fourth-order Runge-Kutta scheme with time res¬ 
olution At = 1.0 x 10 -4 s. The computing box was 
0.1x0.1x0.1 cm 3 . Periodic boundary condition was ap¬ 
plied in the x and 2 directions, whereas solid boundary 
condition was applied at the walls. We reconnected two 
vortices artificially when the vortices approached each 
other more closely than A£. We eliminated vortices that 
were shorter than 5 x A£ = 2.4 x 10 -3 cm. The tempera¬ 
tures were T = 1.9 K (a = 0.21, a' = 0.009), 1.6 K (a = 
0.098, a' = 0.016), and 1.3 K (c* = 0.036, a' = 0.014) 
21]. The initial state consisted of eight randomly ori¬ 
ented vortex rings of radius 0.023 cm. The simulation 
was finalized at t = 2.0 x 10 s for T = 1.9 K and 1.6 K, 
and t = 6.0 x 10 s for T = 1.3 K. 

The quantized vortex tangle in the pure normal flow 
develops to a statistically steady state. Figure[T](b) shows 
the time evolution of the vortex line density L = J c d£/f2, 
where Q is the whole volume. At all temperatures, we 
applied the normal fluid flow with v n = 0.9 cm/s to the 
system, where v n is the value of v n averaged over the 
channel cross section. The value of L develops from the 
initial state and then fluctuates around a constant value. 
This means that QT is in a statistically steady state. 
The statistically steady value of L increases with T be¬ 
cause of a stronger mutual friction. An inhomogeneous 
vortex tangle appears in this statistically steady state. 
Figure [2] shows the configurations of the quantized vor¬ 
tices in the statistically steady state 22], The vortices 
concentrate near the channel walls. A similar structure 
was already reported in the simulations of a nonuniform 


counterflow 17-20J. With increasing T, the inhomoge¬ 
neous structure becomes more apparent and anisotropic 
because of the stronger mutual friction. The configura¬ 
tion and dynamics of the tangle were studied by focusing 
on coarse-grained physical quantities in the latter part of 
this paper. 

The log-law of the superfluid flow is found in the sta¬ 
tistically steady states. The flow direction component 
v*/vo of Eq.([2]) is plotted in Fig. [3](a), where vq is the 
value of Vg at the center y/D = 1 of the channel. The 
values are averaged in the x and z directions and over 
the statistically steady states. The data are plotted as 
a function of the distance y/D from the wall. We can 
observe the logarithmic velocity profile: 


Hi) 


( 6 ) 


where u* is a characteristic velocity of QT, K q is a Kar- 
man constant for QT, and c is a constant value. By fitting 
line to the data, we obtained the values of x*/ K q and c 
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FIG. 2. Snapshots of the vortex tangles in the statistically 
steady states. The left (right) column shows the streamwise 
(side) views of the tangles. 


in Table |T| The quantum Karman constant n q cannot be 
determined. To obtain re q , we should know the physical 
meaning of u*, but we have not constructed a theory of 
the log-law in QT giving u*. 

TABLE I. Numerical results of superfluid velocity Vo at the 
center of the channel, ratio v q / Kq, and parameter c under 
pure normal flow. 


T 

-U0 

Vq/Kq 

c 

(K) 

(s/cm) 

(s/cm) 

- 

1.9 

0.184 

0.141 

1.46 

1.6 

0.079 

0.070 

1.40 

1.3 

0.025 

0.028 

1.14 


Verifying exactly whether the log-law appears may be 
difficult. In fact, researchers of classical turbulence have 
made many discussions of the log-law; for example, some 
researchers argue that the mean velocity obeys a power- 
law 23, 24]. To analyze the log-law more intensively, we 



FIG. 3. (Color online) (a) Flow direction component v a of 
the superfluid velocity v a as a function of the distance y/D 
from the wall. The channel wall corresponds to y/D = 0, 
whereas the center of the channel corresponds to y/D = 1. 
The log-law v]/ = a* [log (y/D) + c] //t q is observed. The solid 
line is fitting line for T = 1.9 K. (b) Another form VKq/Vq = 1 
of the log-law in Eq. is checked. We use the values of 
Vq/nq in Table HI 


check for another form of the log-law: 


_ v q 
V dy K q ’ 


(7) 


which is obtained by differentiating Eq. © with respect 
to y. Figure El[b ) plots Tn q /v q as a function of y/D , 
where T = y(dv*/dy) is the left-hand side of Eq. ©. The 
y-independent regions correspond to the log-law region 
and can be observed in the case of T = 1.9 K and 1.6 K. 
The log-law regions are around 0.2 < y/D < 0.6. The 
log-law region is not wide; therefore, we plan to perform 
another simulation in the future to expand the log-law 
region. As for T = 1.3 K, the y-independent region is 
too narrow to definitively find the log-law. Considering 
the superfluid velocity and the vortex configuration, we 
can classify the whole volume into three regions, namely 
the log-law region Rl, the high curvature vortex region 
Rs near the wall, and the low vortex line density region 
Rc in the center of the channel, as shown in Fig. [2] 

The configuration and dynamics of the vortex tangle 
are important for understanding the log-law in QT be¬ 
cause the superfluid velocity is determined only by the 
quantized vortices. In the following, we report on our 
investigation of the configuration and dynamics of the 
vortex tangle characteristic of the log-law. We define the 
coarse-grained value of a physical quantity Q(£) as 


[QIcg ( r ) — 


u(r)l(r) J c , (r 


Qm, 


( 8 ) 


where w(r) is the local subvolume at r, C'(r) repre¬ 
sents the vortex line within w(r), and l(r) is the lo¬ 
cal vortex line density at r. We divide the compu¬ 
tational box into the subvolumes by using a uniform 
32 x 32 x 32 Cartesian mesh. The separation of the mesh 
is Aa; = 1/32 mm = 3.125 x ICC 3 cm; the mesh satisfies 
Ax > A£. In the following, to focus on the statistical 
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value, physical quantities are always averaged in the x 
and z directions and over the statistically steady states. 



y/D y/D 


FIG. 4. (Color online) Local vortex line density l and 
anisotropic parameter I z as a function of the distance y/D 
from the wall. 

First, the configuration of the tangle was considered. 
The local vortex line density l is plotted in Fig. [4|a). 
The value of l in Rc is smaller than those in Rs and Rl, 
because in Rc the stronger relative velocity between the 
two fluids makes the mutual friction terms in Eq. 
dominant so that the vortices expand rapidly and escape 
from Rc- When approaching the wall from the center, 
the value of l once increases and then decreases near the 
wall. Comparing the data at T = 1.9 K and 1.6 K with 
those at T = 1.3 K, we find that the characteristic be¬ 
havior of the log-law may be the strong inhomogeneity 
of l. On the basis of this observation, the ambiguity of 
the log-law at T = 1.3 K may be attributable to the 
weak inhomogeneity. Figure ®b) shows the anisotropic 
parameter f = [s' • rj] CG , where fi is a unit vector in the 
i direction. A large value of \Ii\ means that the vortices 
tend to be parallel to the i direction. When approaching 
the wall, the value of \I Z \ increases in Rc and Rl and de¬ 
creases in Rs. Although a clear difference is not observed 
between Rc and Rl, we can infer that the log-law needs 
a larger anisotropy than other regions. 



y/D y/D 


FIG. 5. (Color online) Drift velocity V/f and V// as a function 
of the distance y/D from the wall. 

Second, we investigated the dynamics of the tangle 
[2^ . The drift velocity V/j = [s • n] C G °f the vortices 
is plotted in Fig. [5] The data of Vf show that in Rc and 
Rl the vortices move in the flow direction of the nor¬ 


mal fluid component, whereas in Rs the vortices move 
opposite to this direction. In Rc and Rl, the superfluid 
velocity v s is along the normal fluid flow as shown in 
Fig. [3l and in this system, the term including a' in Eq. 
® tends to carry the vortices in this direction. In Rs, 
the value of |u n — v s \ becomes small, so that the mu¬ 
tual friction terms in Eq. ® become less dominant and 
these vortices move against the normal fluid flow by their 
self-induced velocity. A clear difference in Vf cannot be 
observed between Rl and Rc, but between Rs and Rl 
there is difference in the direction of Vf. The data of Vf 
show that the vortices are carried from Rc to Rs by the 
mutual friction. This result and the strong anisotropy I z 
in Rl can cause the momentum transfer of the superfluid 
component to the wall by the vortices. According to the 
analogy with the theory of the log-law in classical tur¬ 
bulence, this mechanism may sustain the log-law in QT. 
When approaching the wall, in Rc the value of || in¬ 
creases, and in Rl and Rs the value decreases. Contrary 
to Vf, a clear difference in the first derivative of Vf is 
observed between Rl and Rc- These results imply that 
both Vf and Vf are important in the log-law region. 

In summary, we demonstrated numerically the loga¬ 
rithmic velocity profile of superfluid flow in QT under 
pure normal flow in a channel. To understand the log- 
law from the configuration and dynamics of the vortex 
tangle, we investigated the coarse-grained quantities of 
the vortex line density, anisotropy, and drift velocity. In 
the future, we will construct a theory of the log-law in 
QT and perform another simulation to obtain a wider 
log-law region. 
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